Exponential estimates of slow manifolds 



K. Uldall Kristiansen and C. Wulff 

Department of Mathematics, 
Technical University of Denmark, 
2800 Kgs. Lyngby, 
DK 
and 

Department of Mathematics, 
University of Surrey, 
Guildford, GU2 7XH 
UK 

Abstract 

In this paper we prove exponential estimates of slow manifolds in analytic systems. 
The results are obtained for general slow-fast systems and finite-dimensional Hamilto- 
nian systems. For general systems we consider the motion in a Banach space with an 
unbounded fast vector-field, while for Hamiltonian systems we consider finitely many 
fast and slow variables. We will prove some conjectures of MacKay from a 2004 refer- 
ence, and the methods we use are based upon the ideas presented in this paper. The 
method does not notice resonances, and therefore we do not pose any restrictions on 
the motion normal to the slow manifold other than it being fast and analytic. 
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1 Introduction 



Singularly perturbed systems involving different time and/or space scales arise in a wide 
variety of scientific problems. Important examples include: meteorology and short-term 
weather forecasting [15j [El |26] , molecular physics and the Born-Oppenheimer approximation 
[20] . chemical enzyme kinetics and the Mechaelis-Menten mechanism [21], predator-prey 
and react ion- diffusion models [22], the evolution and stability of the solar system [131 E3] 
and in the modeling of tethered satellites [271 121] • These systems can also be "artificially 
constructed" by a partial scaling of variables near a bifurcation [25]. The main advantage of 
identifying slow and fast variables is dimension reduction by which all the fast variables are 
"slaved" to the slow ones through the slow manifold. Dimension reduction is one of the main 
aims and tools for a dynamicist and the elimination of fast variables is very useful in for 
example numerical computations. In fact, since fast variables require more computational 
effort and evaluations, this reduction often bridges the gap between tractable and intractable 
computations. An example of this is the long time (Gyears) integration of the solar system, 
see [TU IH]. See also j3j [7J E] for a numerical treatment of slow-fast systems. 

A conjecture by MacKay. 

In [T7J R. S. MacKay conjectured on exponential estimates of slow manifolds in analytic 
systems. It was sketched how slow manifolds with non-degenerate normal behavior could be 
improved and it was then claimed that an iteration would lead to slow manifolds where the 
vector-field would have an exponentially small angle to the tangent-space. To explain the 
method and introduce slow-fast systems let us consider the system 

d t w = eW(w, z), d t z = Z(w,z), (1.1) 

with e a small parameter. The analytic vector-fields W and Z will in general also depend 
upon e, but we shall suppress this dependency throughout. We assume that Z(w, z) = has 
a solution z = ((w), and introduce zq through z = ( + zo to transform these equations into 

d t w = eWo(w,z ), d t z = Z (w,z ) =p {w) + A (w)z + R (w,z ), (1.2) 

with 

Po = -ed w (W(w,()i 
A (w) = d z Z(w,()-ed w (d z W(w,0, 

Ro(w, z ) = / (1 - s) (d 2 z Z(w, ( + sz ) - ed w (d 2 z W(w, ( + sz )) z 2 ds, 
Jo 

and VFo^-So) — W(w,( + ^o)- Here d w is used to denote the (Frechet) partial derivatives 
and we will continue to use this symbol regardless of what object is being differentiated. 
In slow-fast systems one often connects (11.21) with the system 

d T w = W (w,z ), ed T z = p (w) + A (w)y + R (w, z ), (1.3) 
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related to (11.20 by the scaling 



t = et. (1.4) 

If we formally set e = in (11.20 and (ll.3p . then two limit systems are obtained. In (II. 3p . 
the formal limit is singular leading to the algebraic equation: A (w)z + R (w, z ) = (since 
Po = 0(e)). Due to the singular nature of this limit one also often refers to the theory as 
singular perturbation theory. The set of points M = {zq = 0} satisfying these equations is 
called a slow manifold (or critical manifold) and the corresponding system: d T x = Wq(w,0) 
is called the slow subsystem. On the other hand, in ( 11.21) . the formal limit leads to the fast 
sub-, or frozen, system: d t z = A (w)z + Ro(w, z ) with x now considered as a parameter. 
In this system, M is a set of equilibria and the linearization d T 5z = A (w)Szq about these 
determine the classification of the slow manifold. In particular, if zq = is an elliptic or 
hyperbolic equilibrium, then the slow manifold M is said to be normally elliptic respectively 
hyperbolic at the point (w,0). 

One of the main tasks in singular perturbation theory is to determine the fate of the slow 
manifold Mq for e > but small, and thus connect the apparent two different limit systems 
for f / 0. When Mq is hyperbolic then there exists a perturbed, invariant slow manifold M(e) 
for e 7^ nearby [U |5]. For general non-Hamiltonian systems hyperbolicity is generic. On 
the other hand, normally elliptic slow manifolds, which are generic in Hamiltonian systems, 
are unlikely to persist because typical perturbations are believed to destroy them [T7J [18] . 
For these types of slow manifolds, one therefore usually aims for something less: almost 
invariance. For singular perturbed Hamiltonian systems with only one fast degree of freedom 
and an analytic Hamiltonian function, Gelfreich and Lerman [8] showed, using an averaging 
procedure, the existence of an almost invariant slow manifold nearby. By "almost" it is 
understood that the error field (the normal component of the vector field restricted to the 
slow manifold, [17]) is of order 0(e~ c ^ e ), C > 0. For the Hamiltonian example 

H=^x 2 + ^y 2 + v + eyf{u), (1.5) 

with f{u) = Y^=i e ~ n sin(rax) and u = dx /\dy + e~ l du Adv, Neishtadt showed that the slow 
manifold cannot be improved beyond such an estimate, see [8]. The exponential estimate is 
therefore the best one can aim for in a general setting for normally elliptic slow manifolds. 

The method of averaging used in [8] does not extend to several fast variables primarily 
due to the general lack of control of resonances between the fast variables. 

The averaging method of [8] aims at more then we do: the results of [8] do not only provide 
exponential estimates of a slow manifold, they also provide a (9(l)-foliation, parametrized by 
the action variable, of almost invariant slow manifolds. The method therefore, in some sense, 
also addresses stability, not only existence of the slow manifold. The reference [19] extends 
the results of [8] to infinite dimensional slow dynamics. The results of [19] hold true for 
spatially Gevrey smooth solutions, which allow for a Galerkin approximation that separates 
the vector-field into a bounded one and an exponential small remainder. Matching up the 
error from the averaging procedure with the error from the Galerkin approximation the 

authors obtained a slow manifold with O ^exp {—Ce~^\\ error field, where p is a positive 
parameter depending on the Gevrey space. Their results hold true for both Hamiltonian 
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and general systems. The references [29j [30] also provide exponential estimates of particular 
slow manifolds in geophysical models by obtaining optimal truncations of the "super-balance 
equation" (invariance equation) of Lorenz [16]. An extension to more general systems is, 
however, not provided. 

MacKay's method for improving a slow manifold. 

The method of MacKay does not separate normally hyperbolic from normally elliptic slow 
manifolds. The general assumption is just that 1 1 C'^ 77 ) 1 1 1 * s bounded so that the normal 
motion is truly fast. Setting d t z = in (jl.2p therefore gives, by applying the implicit 
function theorem (bearing in mind that i? is quadratic in z ), a solution z = (o(w) close 
to Aq(w)~ 1 po(w) . The graph zq = (o(w) will be the improved slow manifold. To show that 
this is indeed an improved slow manifold, one then straightens out the new slow manifold 
by introducing z\ through zq = z\ + Co- Then the equations become 

d t w = eWx(w, zi), d t zi = Z x {w, z±) =pi(w) + A 1 (w)z 1 + Ri(w, Zi), 

with 

pi = -ed w ((w)W {w,(o), (1-6) 

and so, as Co = 0(e), we have pi = 0(e 2 ) which is the measure of the error field, an im- 
provement from 0(e) to 0(e 2 ), so M\ = {z\ = 0} is an improved slow manifold. MacKay's 
method, though viewed slightly differently, is actually identical to the method suggested by 
Fraser [Sj. An asymptotic analysis of this method is also given in [T2], they do, however, 
only show formal estimates and do not obtain exponential estimates. Furthermore, they do 
not consider the Hamiltonian case. We note that p\ actually vanishes at a true equilibrium 
where Wo(wXo( w )) — 0, and the improved slow manifold Mi = {z\ = 0} therefore includes 
all nearby equilibria. One can continuously perform this procedure as the order of differen- 
tiability of X and Y allows, obtaining a slow manifold with an 0(e n ) error field, but this 
estimate is not in general uniform in n. For analytic systems, however, MacKay conjectured 
that one could obtain an 0{e~ c ^ e ) estimate by applying Neishtadt-type estimates [23l [24] . 
which have also successfully been applied in pij Ej. This is the first main result, which we 
present formally in Theorem 13.11 and prove in section [H In fact, MacKay believed one could 
obtain a stronger result. He conjectured, we guess based on the expression for p\ in (11. 5J) , 
that the error could be estimated point wise in x at each step so that formally an estimate 
of the form 0((e\\W n (w, 0)f) could be obtained leading to an estimate C(e- c/(e||iy(,u ^ (uj)ll) ) 
with n = 0(e^ 1 ) where z = C*( w ) is th e improved slow manifold. We will show that this is 
wrong by considering the following counter example: 

Example 1.1 We consider the simple two-dimensional example: 

w = eW(w, z) =ef(w), z = Z(w, z) =ew — z. (1.7) 

Here z = is actually normally hyperbolic and there is a true invariant manifold nearby 
but this is irrelevant for the purpose of illustrating why MacKay's conjecture is in general 
incorrect. The system has an equilibrium, if f(w e ) = 0, given by (w e , ew e ). We shall assume 
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that the root of / is simple so that f'(w e ) 7^ 0. Note that (jl.7j) is already of the form 
( II. 3p . i.e., z = Zq, Z = Z . The first step of the procedure is then performed by solving 
Z (w, zo) = with respect to zq giving z = £o(w) = ew. Therefore z\ = z — £ solves 

z x = Zi(w, zi) =pi - zi, 

with 

p 1 = Z 1 (w,0) = -e 2 f(w). 
But at the next step when solving Zi(w, z{) = for z\ = (i{w) we obtain 

CiM = Pi = -e 2 /(^), 

and so 

Z2 = ii-(i = Pi- (z 2 + Pi) + e 3 f(w)f(w) = e 3 f(w)f(w) - z 2 . 

Notice that p 2 = e 3 f'(w)f(w) cannot be bounded from above by an expression with \f(w)\ 2 
as a factor. This is due to the fact that w e is only a single root of the dominant order e 3 
term of p 2 - Applying the procedure n times will in general introduce /( n_1 ) in the expression 
for p n . The corresponding term will have a factor of /. This is basically just a result of the 
Leibnitz rule from the continuing differentiation of the £ n 's. □ 

We remark that our results on exponential estimates (Theorems 13 . II and 13.21 below) also hold 
true for normally hyperbolic slow manifolds. Since Fenichel's theorems actually guarantee the 
existence of a slow manifolds in this case these estimates may seem somewhat unsatisfactory. 
However, convergence would imply that the slow manifold would be analytic, which is not 
always the case [2]. 

Hamiltonian systems. 

Normally elliptic slow manifolds are believed to be particular interesting in Hamiltonian 
systems as stability here is associated with oscillatory normal behavior. In Hamiltonian 
systems there are invariant manifolds that are not normally hyperbolic. The most obvious 
examples are regular energy levels and KAM-tori in near-integrable Hamiltonian systems. 
These invariant manifolds are however non-symplectic. Here we are interested in symplectic 
slow manifolds on which we can define a "slow" Hamiltonian system. 

MacKay, still in [17] , suggested a separate method for improving slow manifolds in Hamil- 
tonian systems. The proposed method was described as follows: Consider a Hamiltonian 
H = H(p) with symplectic form uj and a slow manifold Mq. Do the following: 

• Compute an orthogonal symplectic foliation F p so that for every p e M and 

u(z, w) = 0, z E T p F p , w E T p M . 

• Let H p — H\p and solve this for a nearby critical point pi = p\(p). 
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• Put Mi = [piip]} and w Ml = w| Ml - 

Then {H\m 1 ^m 1 ) is an improved slow system. However, we believe that this method has 
some drawbacks. First of all, the method requires the computation of a new slow symplectic 
form at each step. In fact, we believe that the reason for suggesting an alternative to the 
general approach in the first place, is that one wishes to introduce transformations that pre- 
serve the symplectic structure. Moreover, MacKay's method also requires the computation 
of orthogonal symplectic foliations at each step. 

We will therefore suggest an alternative method that circumvent these issues. Our method 
is basically a symplectic extension of the general approach outlined above. We will straighten 
out the improved manifold at each step, ensuring that the transformation involved in this 
procedure is symplectic. Here we will make use of a Lemma 2 in [8]. The slow symplectic 
form uj\ m will therefore remain constant throughout the iteration. The result we obtain is a 
symplectic slow manifold with exponentially small error field containing an initially nearby 
equilibrium, as was also conjectured by MacKay. This is the second main result of the 
paper. We will present this formally in Theorem 13.21 and prove it in section |5j In section [6] 
we present a dynamical consequence of our result on the persiststence of a slow manifold 
with exponentially small gaps. 

2 Notation, assumptions and preliminaries 

Let (W, || ■ || w) an d (Z, \\ ■ \\z) be real Banach spaces and Wc = W©iW resp. Zc = Z®iZ 
their complexifications with norms ||u>i + iw 2 ||w c = ll^illw + 1 1 ^2 1 1 w and \\z\ + iz2\\z c = 
II z i II z + 1 1 ^2 1 1 .z- We will from now on denote all norms, including operator norms, by || ■ ||. 
Hopefully it will be clear from the context what norm is used. Then / : Vc — > Zc, with 
Vc an open subset of Wc, is analytic if it is continuously differentiable, i.e., if there exists 
a continuous derivative d w f : Vc —> £(Wc, Z c ), where £(Wc,-2c) is the Banach space of 
complex linear operators from Wc to Zc equipped with the operator norm, satisfying the 
following condition 

\\f( w + h)-f(w)-d w f(w)(h)\\ = 0(\\h\\ 2 )- 

By real analytic we will mean analytic and real when the arguments are real. The higher 
order derivatives can be defined inductively and d™f becomes a map 

dlf:Vc^£ n (W c ,Z c ), 

from Vc into the Banach space £ n (Wc, Zc) of all bounded, n-linear maps from Wc x ■ ■ • x Wc 
(n times) into Zc- See [HI Appendix A] for a reference on analytic function theory in Banach 
spaces. When V is an open subset of W and x > then, as in (8], we define V + iv to be 
the open complex ^-neighborhood of V: 

V + iu = {w E W c | d Wc (w,V) < u}, 

where <i>v c is the metric induced from the Banach norm || • ||. In the following let Bf(z) C 
Zc = {m£ Zc, \\u — z\\ < r} denote a ^c-open ball of radius r > around z in the Banach 
space Zc- We frequently need the following Cauchy estimate: 
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Lemma 2.1 Assume that f : Vc — > Zc is analytic and that f is bounded on B u (wq) C Vc 
for some v > 0. Then 

\\d w f(w )\\< SUPw ^ W0) ll/H " . (2.1) 



Remark 2.1 Let / : V + iv —> 2c be analytic and bounded and v > uq > 0. Then we can 
apply this estimate to any Wo E V + z(^ — ^o) to obtain: 

sup ||^/(w )|| < , 

which we will write compactly as 



l>0 



C.(T 



This is the form of Cauchy's estimate that we will be using. Similarly, we will by 
denote the sup-norm taking over the domain (V + iv) x [S + %o~ ) . □ 

We will also use the following generalized version of Taylor's theorem 

Lemma 2.2 If f : Vc — > Zc is n times continuously differentiable, n >1, and if the segment 
w + sh, < s < 1, is contained in Vc, then 

f(w + h) = f(w) + d w f(w)(h) + ■■■ + —t—d^fMh"- 1 

[n — 1)1 

"1 /1 n \n-l 



l (n-1)! 9U{W + Sh){K " ' ' k)dS - (2 - 2) 



The integral remainder is bounded by ^j- sup 0<s<1 \\d™f(w + sh) \\ . □ 
Here, we write for m E N, 

d™f(w)h m = d™f(w)(h,...,h). 
Assumptions on general system. 

By straightening out the slow manifold z = ({w), the general slow-fast system (11. ip was 
transformed into ( 11. 2ft . repeated here for convenience, 

d t w = eW (w, Zo), d t z = p (w) + A (w)z + Ro(w, z ), 

and we will now state our assumptions on this system: 

(Gl) Let V C W be open and let u > 0. We assume that Ao(-)" 1 : (V + iv ) ->■ C(Z C , Z c ) 
is real analytic with ||A (u»)~ 1 || i/0 < 
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(G2) Let ES C Z be open and bounded, with C z > such that ||zo|| < C z for all zo E 
S+icjQ and let ctq > 0. Assume that po : (V+wo) — > -2c W^o : (V+Wo) x (<S-Hcr ) — >■ Wfc 
are real analytic. Furthermore, we assume that the norm of Wo is uniformly bounded 
by C Wq on (V + iu ) x (S + «cr ), and 

«y = ||po|k = 0(e), (2-3) 

where e > is small; 

(G3) Rq(-, •) : (V + iz^o) x (<S + icr ) — is real analytic and quadratic in z so that 
Rq(w,zq) = '52iJ2j>2 a ij w '' z 0' with proper interpretation of as multi-linear oper- 
ators. Furthermore the norm of Rq is assumed to be uniformly bounded by Cr on 
(V + iv ) x (S + ia ); 

We now write eW as 

eW = (C Wo e)(C^ o W ), (2.4) 

and let CV e be our new small parameter, which we continue to denote by e, and replace Wo 
by C^ q Wq. We will also continue to denote the new Wo by the same symbol, hoping that 
this will not cause unnecessary confusion. We therefore continue with 

C Wo = 1. (2.5) 
Assumptions on Hamiltonian system. 

In section [3721 we will consider a real analytic Hamiltonian system H = H(w,z), where 
w = (u, v)e W, dim W = 2dw, are the slow variables and z = (x, y)E Z, dimZ = 2dz, are 
the fast variables, with d = dw + dz degrees of freedom, <i w slow- and dz fast ones, and the 
symplectic form 

ijj = dx A dy + e~ l du A dv (2.6) 

(see also Neisthadt's Example (jl.5p ). As for the general case above we consider a slow 
manifold Mo set of constrained equilibria: 

M = {z = CM = {(*(u,v),(»(u,v))\d t H(w,z) = 0}. 

This will be our (first) slow manifold which we aim to improve upon. We will straighten 
z = C(w) out by introducing new coordinates through the generating function: 

G (u, v , x, y ) = (x, y ) + e~ x (u, v ) + g (u, v , x, y ), 

with 

g (u,v ,x,y ) = -(C(u,v ),yo} + (( y (u, v ), x). (2.7) 

Here (a, b) = a±bi + . . . + a n b n for every pair a = (ax, . . . , a n ), b = (b\, . . . , b n ) E C n . Then 
(u,v,x,y) (u ,vo,xo,yo), for e sufficiently small, is a symplectic transformation given 
implicitly by the equations: 

x = d yo G = x-C, V = d x G = y + ( y , 
uo = ed VQ G = u + ed Vo g, v = ed u G = v + ed u g. 
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This transforms H into 

1 

H (w , zq) = H(w, z) = h (w ) + (po[w ), z ) + -(A (w)z , z ) + r (w , z ), (2.8) 

with p = C( e ) on (V + ivq) x (iS> + icr ) where V, S are as before and is , cr > 0. We now 
assume the following: 

(HI) A = dz o H \ Zo=0 : S + iao is a symmetric matrix valued real analytic 

function satisfying HAjIUo — Ca and || A (w)^ 1 1|„ < 

(H2) p = 9 ZQ H \ Z0=0 : V + iu — > C(Z, C) and h = H \ zo=0 : V + iuQ ^ C are real analytic 
functions and 5o = \\po\\u = 0(e) is small. We also assume that ||/io|U < Ch and 
||<9„,/zo|L < C' hQ ; 

(H3) The function r : (V + «z^ ) x (<S + 2cr ) — > C, is real analytic, satisfies r = C(zq), and 
its norm is uniformly bounded by C ro ; 

(H4) There exists a locally unique equilibrium (w e , C{w e )) in the (w, z)-coordinates which in 
the transformed (w , 2 )-coordinates takes the form (u>q, 0) G V x S. 

Note that we can assume Zq = in (H4) without loss of generality. Indeed, the equilibrium 
(u> e , £(u> e )) transforms to (wq,Zq) with 

«o = C(w e )-C(«",«o)» 

and so by introducing the affine symplectic transformation (iuq, ^o) ^ (^o> — ^o) we move 
this equilibrium to (wq,0). It is important to start our iteration from Zq = - we can 
then ensure that our improved slow manifolds, which we will define iteratively, contain this 
equilibrium. Obviously we could also transform Wq = but this will not be necessary. 

3 Main results 

3.1 General slow- fast system 

Our main result for general slow-fast systems is the following: 

Theorem 3.1 Consider U.2\) and assume that assumptions (Gl), (G2) and (G3) hold true 
and let v_, q_, £o be chosen positive constants satisfying 

v < V\ < Vq and a < a\ < a where V\ — r/ — £o> °~\ = °o — £o- (3-1) 

Then there exists an e > so that for every e < e the following holds true: There exist 

(i) a transformation (V + iv) x (S + iq) — > (V + iuq) x (5 + iao), 

(w, z*) i y (w = w,z = z* + C(w)), 

transforming $1.2]) into the normal form 

d t w = eW*(w, z*), 

dtz* = p*(w) + A*(w)z* + R*(w, z*), 
where W*(w, z*) = W (w, z ); 



9 



(ii) constants C\, . . . , C5 only depending upon the previous constants so that 

||p*HII < C 1 e 2 ||^( W ,0)||e- [mln2/(4Kie)1 , m = min{^ - v, a x - a}, (3.2) 

where K x = K (l - K qT 2 5 (^a + 2C Rq K { X\ , and 

\\z* — 2o||i/,CT < Ci^-i ||A* — ^olli^ < 
II #* - ^o|| K , 2 < C^, || W* - W ||^ < C 5 e. 

In other words: U.2\) has a slow manifold given by z* = which is invariant up to the 
exponential small error p*, which vanishes at equilibria where W*(w,0) = 0. □ 

3.2 Hamiltonian slow-fast system 

In the Hamiltonian case it is a bit more difficult to be explicit about the dependency of 
new constants on the old constants without introducing an overwhelming amount of clutter. 
Theorem 13.21 will therefore appear slightly less explicit than Theorem 13.11 

Theorem 3.2 Consider the Hamiltonian H = H(wq, Zq) from \2. 8\) and assume that as- 
sumptions (HI), (H2) and (H3) hold true. Furthermore, let v_, q_, £0 be chosen positive 
constants satisfying y_ < v\ < uq and a < o~\ < o~q where V\ = uq — £0 and o\ = o~q — £o- Then 
there exists ane>0 so that for every e < I the following holds true: There exist 

(i) a transformation (V + iv) x [S + ia) — > (V + iv ) x (S + ia ), (w*,z*) ^ (w , z ) 
transforming \2. 8\) into the normal form 

H^w^z*) = H(w ,z ) = h*(w*) + (p^w*)^*) + -(A*(w*)z*,z*) + r*(z*, u>*). (3.3) 

(ii) constants C±, . . . , C5 only depending upon the previous constants so that 

and 

Ik* - ZoL<7, Ik* - w \\ Et2 _ < C 3 e, 
II A - A \\u,a < C 4 e, ||r* - r \\ Et a: < C 5 e. 

(Hi) If (H4) holds, so that (wq, Zq = 0) is an equilibrium, then p*|( w e,o) = 0. 

In other words: 1(2.8) has a slow manifold given by z* = which is invariant up to the 
exponential small error p*, and by (Hi) contains the nearby equilibrium (wq,Zq = 0). □ 

We can now directly address stability: 
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Corollary 3.1 Consider the Hamiltonian H* = H*(w,z) from 113. 3\) on the real domain 
V X S. If Aq is positive definite then L(w,z) = \(A*{w)z, z) +r*(w,z) is an approximate 
Lyapunov function for \\z\\ and e sufficiently small, and there exist constants c\ and C2 so 
that 

\\z(t)\\ < 0(e~ Cl/e ) for 0<t< c 2 e" 2 , 

when z(0) = 0. 

Proof By assumption A* is positive definite for e sufficiently small, and hence there exist 
constants A> and A> so that 

< Alkll 2 <L< A||z|| 2 (3.4) 

for ||z|| small and w G V + iv. Let J z , J w be such that the symplectic form u from (12.61) 
satisfies 

u((z 1 ,w 1 ),(z 2 ,w 2 )) = (zi, J~ 1 z 2 ) +e~ 1 {wi, J~ 1 w 2 ). (3.5) 
When differentiating L we then obtain, 

L = (d z L, J z (p* + A*z + d z r*)) + e(d w L, J w d w H) 

= (A*z + d z r*, J z (p* + A*z + d z r*)) + e(d w L, J w d w H) 
= (A*z + <9 z r*, J z p*) + e(d w L, J w d w h*(w*)) 
<C 6 e 2 e- C2/t + C 7 e sup L(w,z(t)) : 

w&V+iu 

on (w, z) G V x S for some constants Cq and C 7 . Here we have used a Cauchy estimate on 
d w L. Integrating this inequality from t — to t and using (13.41) we find that 

X\\z\\ 2 < L < C 6 e 2 te~ C2/t + C 7 e [ sup L(w, z{s))ds < C 6 e 2 te~ C2/t + C 7 Xe [ \\z(s)\\ 2 ds. 

JO w Jo 

We have here used that L(0) = since z(0) = 0. Then by Gronwall's inequality in integral 
form [1] we obtain 

Ikll 2 < CeX-'eHe-^e ^ 1 *, 
and therefore while < t < C 2 \/(2C 7 \e 2 ): 

\\z\\ < 

" " ~ V 2CVA 

completing the proof. ■ 
Note that this upper estimate 0(e~ 2 ) is large, even on the fast time scale r (II. 4p where it is 
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4 Proof of Theorem 13.1 



We will first state and prove two key lemmata. The first lemma will be an application of the 
implicit function theorem, while the second one will be "The Iterative Lemma". We will use 
these lemmata successively, so assume that p, A and R (in place of po, Aq and Rq) satisfy 
the assumptions (Gl), (G2), and (G3). 

Lemma 4.1 Suppose that < k < a — Kb where 

8=\\ P {w)\\ v <k 2 K- 2 C^. (4.1) 

Then the equation 

= p(w) + A{w)z + R{w,z), (4.2) 
has a locally unique solution z = ((w) G Z c satisfying: 

HCHll<iq P H||, (4.3) 

for every w G V + iu. Moreover ((w) is analytic in w G V + i 
Proof Re-arranging ( 14.21) and applying the inverse A(w)~ l gives 

( = -A(wr\p(w) + R(w,()). (4.4) 

Put (o(w) = —A{w)~ 1 p{w) and consider ( = ( + z so that 

z = -A(w)- l R(wXo + z). (4.5) 

Note that ||Co|| < f ||p(w)||. We will denote this upper bound by r := y||p(u>)|| and highlight 
r < f 5. Set 

T(z) = -A(w)- 1 R(w, Co + z). 
We will show that T is a contraction on -B^(O) C S + ia. By Taylor's formula if follows that 

R(w,z)= / (1 — s)d 2 z R{w, sz)z 2 ds, 
Jo 

and therefore by applying a Cauchy estimate we obtain for ||z|| < K5 

\\R(w,z)\\<l sup II^Ks^llll^ll^^PII 2 . 

^ 0<s<l K 

Here < k < o — K5. Therefore 



\T(z)\\ <^{2r) 2 <K 2 ^5r<r 

2 K z K z 



using that 2r < Kb and the assumption (14. ip . and hence T : Bf{Qi) — > Bf(0). Next, we use 
Taylor's formula to obtain 



d z R(w,z)= / d 2 z R{w,tz)zdt, 
Jo 
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and then apply a Cauchy estimate to bound d z R: 



9C 

\\d z R(w,Co + z)\\ <— % + z||. 



Therefore 



K 2 ^§5 < 1 



K 

using (14. ip . This shows that T is a contraction on -B^(O) and there exists a unique fix point 
z{w) of T. In particular, C(^) — Co(^) + z(w) solves (14 .4p and HCC 10 )!! < 2r = _?T||p(u>)||. By 
[TUl Section 1.2.6] the map ( : V + iv — > Z is analytic (this is where we need w — > A~ l (w) 
to be analytic). 

Lemma 4.2 (The Iterative Lemma) Let 

£ > 2Kb (4.6a) 

and 

S < ^k 2 K- 2 C r \ e < 2k 2 , (4.6b) 
v — £ > and < « < o — £. T/ien t/ie transformation 

(w,z+) i ^ (w,CO) + 
/rom (V + — £)) X (<S + «(<x — 0) * n ^° + Zl/ ) x (5 + i(cr — £/2)) transforms 

d t w = eW(w, z), 
d t z = p{w) + A{w)z + R(w, z), 

into the normal form 

d t w = eW + (w,z + ), 

d t z + = A + (w)z + + p+(w) + R + (w, z + ). 



Here \\p\\ x < 6, ( is the solution from Lemma 4-L we define W + (w, z + ) = W(w,C + z + ) and 
the following estimates hold true 

llp+HII < y5||W + ( W ,0)||, (4.7a) 
IKH- 1 !! < j(l-^K~ 2 (^K + 2C R K^y\ (4.7b) 
\\A + -A\\ v _t < K - 2 8(^-K + 2C R Ky (4.7c) 

\\R+-R\\ v -te-i < € ^5(2 + K~ 1 C z )+C R ^^ + K- 2 K5(K5 + 2C z )y (4.7d) 

for w G V + i{y — £) and where ||y4(w) _1 ||„ < y , ||-R||^cr < C_r awe? ||z|| < C z . In ■particular, 
= \\P+\\v-£ satisfies 

5 + < ^5. (4.8) 

□ 
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PROOF Applying the transformation z = ( + z + gives 

d t w = eW + (w,z + ), 

d t z + = A + (w)z + + p+(w) + R+(w, z+), 

and using the fact that ((w) solves ( 14. 2ft we obtain 

A+(w) = (A(w) - ed w ((w)d z W(w, CM) + d z R(w, ((w))) , (4.9a) 

P+ (w) = -ed w ((w)W(w, CM), (4.9b) 

W+{w,z + ) = W{w,C + z+), (4.9c) 



and 



R+(w, z+) = -ed w C(w)W(w, z + ) + R(w, z+), (4.10) 
W(w, z+) = W(w, ( + z + )- W(w, C) - d z W(w, ()z + , (4.11) 
R(w, z+) = R{w, ( + z + )- R(w, C) - d z R(w, Qz + . (4.12) 



We have 



e K 

^\\W{wX)\\<-^5\\W + {^Ml (4-13) 



for w G V + i{y — £). Here we have used the Cauchy estimate: 



\\dA-i<^f, (4-14) 

for v — £ > and £ > and Lemma 14.11 To complete one step of the iteration we also have 
to estimate 

A+(w) - A{w) = -ed w (d z W(w, + d z R(w, (), 
and R + — R appropriately. For the former first note that by (I4.14p for w G V + %{y — £) 
\\d w (d z W(w,C)\\c(z,Z) < \\d w (\\c(w,z)\\d z W(wX)\\c(z,w) 



< —K . 



using another Cauchy estimate, noting (12. 5p and cr — £ > k > and ||C||^ < £/2 by Lemma H~T1 
and (I4.6al) . Moreover 



\d z R(w,C)\\ = II f d 2 z R(w,s0s(ds\\ < 2C r k~ 2 
Jo 



for all w G V + i(z/ — £), and therefore for all such w, as ||C||f < by Lemma [4. 11 



\A+ - A\\ u ^ < —5k^ + 2C r k~ 2 K5 = k~ 2 5 1—k + 2C R K ) . 
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Hence, 



< \\A(w)z\\ 



= \\(A + (w) - (A + (w) - A(w))z\\ 

< \\A+(w)z\\ + \\A + {w) - A{w)z\\ 

< \\A + (w)z\\ + \\A + (w)-A(w)\\ u ^\\z\ 

< \\A+{w)z\\ + k~ 2 5 (^-k + 2C R K 



so that 



if 



| (l - (j-k + 2C R K ) ) \\z\\ < \\A + (w)z\ 



KK ~ 2 . 5 ( e Ji K + 2C R K] <1. 



This is true due to (14. 6|) : 

^f-S (^-k + 2C R K) < ^ + K*k-*C r 5 

1 1 , 
< - + - = 1, 

2 2 

We highlight that by (14.6b|) this choice of 6 satisfies the condition (14. ip required in Lemma |4~2 
- it is in fact a factor 2 stronger. Hence 

P+M'l < § (i - ^-pW 1 + 2C f 

for all w G V + i{y — £)■ 

For R + - R, with R + from (I4.10p we first estimate W from (14. lip and # - i? from (I4.12p . 
This gives: 

u,e-t < 2 + K~ l C z 
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and 



\\R - R\\ v -t, ff -t < \\R{w, ( + z + )- R(w, z+) Wv-fr-i + \\R{w, () \\ u -(,<,-t 
+ \\d z R(wX)z+\\v-z,< r -t 



<\\ [ d z R(w,sC + z+)dsC\\ v -t,*-t + \\ [ (l-t)d 2 z R( 
Jo Jo 



WitC)C 2 dt\\ u -z, a -£ 



d 2 z R(w,t()(z + dt\\ u _^ 



< II^IU-^IICL + ~\\C\\l + ^IICII.c 

2 K z K z 

2C 

< —jr-K5 + C r k~ 2 K5 (KS + 2a) 
= C R + k~ 2 KS (KS + 2a) 



Therefore 



\\R+ - R\U-s,*-t < e\\d w C\l^\\W\l^ >(7 ^ + \\R- R\\^s 

< ^-5(2 + k~ x C z ) + C R + k~ 2 K5 (KS + 2C Z 

Following this lemma, we introduce new constants: 

K+ = k(i- (^-k + 2C r k] ) _1 , (4.15a) 



2 v £ 

C R+ = C R (l + ^- + k~ 2 KS (KS + 2cSj + ^p(2 + k" 1 ^), (4.15b) 

(7+ = O- - £, 

v+ = v-i. 

We will now apply Lemma 14.21 to (ll.2p once with £ = £ = C(l) < max(z/ — v, <7 — a) 
lar<?e. Due to ( 12. 3ft and (14.7aj) there is a C5 > such that 

5x=C 5 e 2 . (4.16) 

So we start from (V + iv\) x (S + zoi) with i^i = z/ — £0 and cii = a — £ - By ( 12.3ft and ( 14.3[) 
we have \\zi — z \\ u < K S. Moreover we can bound A^ 1 and Rx by 

2x \\ ~ 1 



IIVIU < *i = JT (l - ^ + 2^0)) 



resp. 



pilkai < = Cifo ( 1 + H' 2 W (Wo + 2C 2 )1 + ^^(2 + a- 1 C z ), 

\ £0 / £0 
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taking ft = a in (14.15aj) and ( 14.15b[) (noting that by (13. ip this choice of k satisfies k < 
a — = fi as required in Lemma [4.21) . We then continue to apply Lemma [4.21 with 

a = 2K n e, (4.17) 

for n < N, where N is such that vn — £,n > U. an d &n — £jv > <Z- Here we have defined 



V n+ 1 = Uy - ^2 6 
1=1 

n 



i=i 



Note that the form of £ n from (I4.17P allows us to apply "The Iterative Lemma 14.21 '. since at 
each step 

£ n > 2K n S 1 = 2K n Cse 2 > 2K n 8 n 

provided e is sufficiently small so that e < Cj 1 and provided that 5 n < 6\. Here we used 
dUB}. Then by (Ojl . 

if n <5„ n AT n x • • • x K x 

< 2-%. (4.18) 

Here <5 n+ i = \\p n +i \\v n+1 - To prove the theorem, it is left to be shown that £ n = 0(e) with an 
order constant which is bounded in n so that by choosing N = 0(e~ l ) iteration steps the 
estimates in Theorem 13.11 apply. We therefore estimate K n+ i and C Rn+1 : 

K n+l <K n (^-^a-H n {^n + 2C Rn K^ * 

<K n (l- 2~ n+1 K n qT 2 5 1 (\a + 2C Rn K n )) , (4.19) 



C Rn+1 < C Rn (l + + ar 2 KJ n (KJ n + 2C Z )) + ^(2 + k' x C z ) 

< C* n (1 + 2-" +1 e- 1 <5 1 + 2-V- 2 if n( J 1 (2""^^ + 2C Z )) + 2"% (2 + a^C*). 

(4.20) 



If 



0-%(^ + 2C^) <— (4.21] 



and CKj = J^- then (I4.19P implies 

«n+i < (1 - 2- n ~ 2 a n )' 1 a n . (4.22) 

It is easy to show that 



n+l 



a n < J— < 2, (4.23) 
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and so 

K n < 2Kt for all n < N. (4.24) 

This estimate could be improved, but we do not need to do so. Similarly for C Rn , we set 
Pn — an d obtain 

0n+i <Pn{l + 2~ n+1 e- 1 5 1 + 2- n+2 qr 2 K l 5 1 (2- n K l 8 1 + C Z ))+ 2~ n S 1 
with 5i = C^%(2 + qT x C z ), using the estimate (I4.24p obtained for K n . Let 

cT 2 K\e{K\b~\ + C z ) < X - (4.25) 

so that 

n+1 < (1 + 2~ n+2 C s e)(3 n + 2~ n h. (4.26) 

Let N G N be such that /3 n < 2 for n < N . Then for n < N we get, using the geometric 
series, that 

n </3 1 + 8C s e + 5 1 . (4.27) 

Using ( I4.13P and choosing e small we see that (3 n < 2 for all n or simply 

C Rn <2C Rl . (4.28) 

Returning to our assumptions, (I4.2ip and (14.251) . we realize that due to (I4.16P it suffices to 
take e so that the inequalities 



2=~ 1,1 7 - 8ATi 
<T*K x e(K x C s <? + C Z )<±, 

are satisfied. Notice also that the hypotheses of Lemma 14.11 are satisfied at each step. 
Therefore, by {Qfl and (04]) . 

< 4Kie. 

and so we can perform 

m 



N 



AK x e 



m = minjo"! — g_, V\ — u}. 



iterations before one of the conditions fjv+i > v, <Jn+i > Q_ is violated. 

If we inspect p^ + \{w) further we realize that the error field estimate can be made point- 
wise in x improving near actual equilibria: by (I4.13P and (I4.18P we have 

\\Pn+i{w)\\ < ^-^6 N \W N (w,0)\ 

< 2- [m/ ^ e)] S 1 \W N {w,0)\ 

< C s e 2 e- [mln2/{iKie)] \W N (w,0)\. 
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Finally, provided e is sufficiently small, there exists constants c±, . . . ,04 depending only upon 
the previous constants so that by (IP]) . (OS) . flOl) and fl4~T6j) 



JV N 

\\z N+ i - Zl \\„< < J2 2 ~ jRiSi ^ cie2 - ( 4 - 29 ) 
3=1 3=1 

Similarly, by (14. 7c p . using (I4.17P and that K n and Cr u are bounded independently from n 
(see (jjjlD and ( EHD ), fl4T8|) . and ( 14161) we get 

H^at+i - < c 2 e 2 . 

Moreover, using (l47d1) . (OTjl . f l4T24j) . f l4T28|) . (OS]) , and fl4TT6|) for the first estimate and the 
mean value theorem and (I4.29P for the second estimate we get 

H-Rjv+i - Ri\\ E ,a < c 3 e 2 , \\W N+1 - Wi\\ K> a: < c 4 e 2 - 
Since \\z\ — ZolU is ^( e ) by Lemma [4. II and (G2) the theorem follows. 

Remark 4.1 Again we highlight (cf. Example 1 1.1 ft that it is in general not possible to obtain 
an estimate of the form 0(e~ c ^ t ^ w *^ w '° s> ^) , since at each step of our iteration we must apply 
the Cauchy estimate (14.141) forcing us to estimate £ on V + iv - through (14.31) this gives rise 
to the factor CV = 1 rather than ||W^(w,0)|| at each step. □ 



5 Proof of Theorem 3.2 



The proof of Theorem 13.21 is again based on an iterative lemma. As we will use this lemma 
successively we will in the following consider the normal form Hamiltonian 

H(w, z) = h(w) + (p(w), z) + -(A(w)z, z) + r(w, z), (5.1) 

(in place of H from (12. 8p ) satisfying the assumptions (HI), (H2), (H3) and (H4) on (w, z) G 
(V + iv) x (S + ia), with 5 = \\p(w)\\ v small. Let z = ((w) be the solution of 

d z H(w, z) = p{w) + A{w)z + d z r(w, z). 

If ||A(w) _1 || < y then by Lemma [4.11 this solution satisfies the estimate \\C\\u < KS for 

5 < \k*K-*C;\ (5.2) 

where the "extra" factor |k compared to above in Lemma [4.11 comes from estimating d z r = 
0(z 2 ) rather than just R. Let C{w) = (( x (w),( y (w)) an d consider the generating function 
G(u, v + , x, y + ) from (12.71) . We then consider the symplectic transformation (w + , z + ) \-> (w, z) 
defined implicitly by the equations 

x + = d y+ G = x - C(u, v+), y = d x G = y+ + ( v (u, v+), 
u + = ed v+ G = u + ed v+ g(u, v + , x, y + ), v = ed u G = v + + ed u g(u, v+, x, y + ). 
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Lemma 5.1 The function g from {2. 7| ) satisfies \\g\\ u ,a < CollClk? where Co only depends 
on the domain (V + iv) x [S + io~). □ 

Proof This is straightforward and follows directly from the estimate 

lkll<||C1ll|y + ll + lieillkll<c D ||ci|. 

We will write this transformation as 

z = z + + a(w + , z + ) =z + + C(^) ; w — w + + e P{ w +i z +), (5.3a) 
with z + = (x + ,y + ), w + = (u + ,v + ), w = (u,v + ), and use 

P = Po(w+) + (Pi(w+),z+) + ^(f3 2 {w+,z + )z + ,z + ), (5.3b) 

and 

w = w + + e0o + 0i, z + ) + -(02Z+, z+)). (5.3c) 

Here ft = (ft, 0) when f3 l = (ft, ft). 
We have 

Lemma 5.2 Let g satisfy \\g\\ u ,cr < and let £ be such that v — £ > and a — £ > 0. 
T/ien (u>+, z+) i— ?• (w;, z) is a symplectic transformation (V + i[y — £)) x (S + i(cr — £)) — y 
(V + i(v - £/2)) x (S + i(a - £/2)) satisfying 

\\a\\ v ^, 11/%,-^ < < | (5.4) 

□ 

Proof This is basically Lemma 2 in [8]. For comparison note that their es, ea, e(3 is our g, 
a resp. 0. ■ 

In particular 

4e|| sll 

e llA)||i/-f = \\ w + - w ( w +,®)\\v-t < — £ — , ( 5 - 5a ) 
and using Cauchy-estimates on c^ + /3(u> + , 0) and c?2 + /3(w + ,0) we get 

\Wi(w + )\\ < K^m^, \Mw + ,o)\\ < 2K- 2 m\^- (5.5b) 

Corollary 5.1 The estimate for ||a|| I/ _^ cr „^ in ( [5.^1 ) can 6e improved to 

\\a\\ u ^ a ^ < \\(\\ u . 
Proof By the previous lemma w\ v -^ jCr ^ C V + i(v — £/2) and so 

II«IIi/-€,«t-{ = IK(*)IU-^-« < ||C||^/2 < 
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The fact that (w+, z+) (w, z) is well-defined with domain (V+i(v — £)) X (S+i(a — £)) and 
co-domain (V + %{v — £/2)) x [S + %{p — £/2)) was crucial here and will be in the following. 
The v — ^jl and cr — £/2 terms in the co-domains allow for a step of £/2 to apply Lemma 12.11 
to estimate derivatives on (V + i(V — £/2)) x (<S-M(cr — £/2)) by function values on the larger 
domain (V + iv) x (5 + io - ), c.f. the Cauchy estimate (12. ip . This introduces a factor of 2£~ 1 . 
We are now ready to state and prove the "Iterative Lemma for Hamiltonian systems": 

Lemma 5.3 (The Iterative Lemma for Hamiltonian systems) Consider 115. 1\) and 
assume that the assumptions (HI), (H2) and (H3) hold true, so that h, A and r satisfy 

\\h\\ v < C h , \\d w h\\ uH/2 < C' h , 

K 

\\A\\i> < C a-, \A ||p ^ "2") 



Furthermore, let 



ML < C r . 



f 2 1 



f > i\~C D <5 < \, 5< -K~ 2 C; L ti\ z/ - £ > and < k < a - £. (5.6a) 

8 3 

Tnen inere exists an e and constants Ch.i, Ch,2, c&, ca,i, c a,2, Cr,x and c r ^, depending only on 
the previous constants, so that for e < e the symplectic transformation \& + : (w + , z + ) 1— > (w, z) 
from Lemma VxJ^ mapping (V + i(v — £)) x (S + i(a — £)) into (V + z(z/ — £/2)) x (S + i(a — £/2)) 
transforms H = H(w,z) from \5. 1\) into 

H + (w + , z + ) = H(w, z) = h + (w + ) + (p+{w + ), z + ) + ^(A + (w + )z + , z + ) + r+(w+, z+), (5.6b) 

\\ h + ~ ti\\u-t < c h ,iO + — ^— , ||9«,ft+||i/-3e/2 <c7 fc + — 1 ^— , Hp+IU-c < 

°a,i^ 11 I, c- , c r , 2 e5 

I 5 



|| A + - A\\ v -e < c A>1 S H ^— , ||r+ -r|| v _ €|ff _£ < c r>1 5 + — ^— . (5.6c) 



In particular, 



\\(A + )-%-t < ^ := f (1 - f + c A2 ) 5)-, (5.6d) 



and <5 + = Hp+ll^t satisfies 



5+ < ^-5. (5.6e) 



□ 



Remark 5.1 The constants c^i, 0^2, c$, ca,i, ca, 2 and c r> i, c rj2 depend polynomially on K, 
Ch, Ch, C' h , Ca, C r and S. This is made more explicit in the proof below. D 
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Proof It is convenient to introduce 

w + = w\ z+=0 =w + + e(3 (w + ), (5.7) 

and 

(s = ((w + + se(p -Po)), (5.8) 

with s G [0, 1], so that in particular ( Q = ((w + ) and £i = ((w + + e A))- Note also, due to 
(jOSjl and (Ell, that 



Ci = CH* + =o). (5.9) 

We will use this frequently in the following where we expand the new Hamiltonian H + (w + , z + ) 
H(w,z) about z + = to put it in the normal form ( 15.6bj) with 

h + (w+) = H+(w+,0), P+ (w + ) = d z+ H + (w + ,0), A+(w + ) = d 2 z+ H + (w +} 0), 

and 

r+ = H+ - (h+ + (p+, z + ) + -(A+z+, z+». 

Due to (15.71) and (15.91) we obtain 

h+(w+) = H + (w + , 0) = H(w + , d) (5.10) 

= H(w + , 0) + / d z H(w + , sd)Cids, 
J o 

= + e/3 ) + / d z H(w + , sCi)Cids, 
Jo 

= h(w + )+ / (e9 w /i(u7+ + se/3 )/3o + 5 2 ^(^+,sCi)Ci)^. 



Using that by (I5.6al) and Lemma |4~T1 we have £ > 2K5 > 2||£ L || i/ _f/ 2 an d employing (I5.5al) 
we obtain 

- h(w + )\\ v -£ < e\\d w h\\ u ^ /2 \\^oh~^ + ||d*#||«/-€/2,e/2||Cl|i/-£/2 
< ^^y^ + k-'ChKS 

<c M d + — ^— , (5.11) 

introducing the constants 

In (15. lip we have used (I5.5ap , Lemma 15.11 and Lemma 14. 11 
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Moreover, using ( 15. 3 aft and fl5.3cj> to compute J^-| z+ =o and -^\z + =o we obtain 
p + (w + ) = d z+ H + (w +: 0) = d z H(w + , Ci)(/ + ed w CJi) + ed w H(w + , £)& 

= e / dlH{w + Xs)d w Uh-^)ds{I + ed w Ch) + ed w H{w + Xi)^ (5.12) 
Jo 

where we have used that d z H(w, ((w)) = to arrive at 

d z H(w + , Ci) = d z H(w + , ((w + + e(p - A,))) = d z H(w + , ((w + )) +e / -d z H(w + , ( s (e))ds 

s v ' Jo de 

=0 

= e / d 2 z H(w + Xs)d w (s0o-Po)ds. 
Jo 

Also by definition 

d w H(w, z) = d w h(w) + (d w p(w), z) + -{d w A(w)z, z) + d w r(w, z), 

and therefore, using that r is cubic in z by (H3), we can estimate the factor of e/3i in the 
last term in ( 15.121) by 

H^tfodH^ < H^^H^/a -h ||^ H ,/»|| 3e _e/2||Cll 3 ^-€/2 -I- Ill^^C^llxr-^^illClS-^/a 
+ ^ll^|| x - 5A5 / 2 ||Cll^/ 2 



<^ + 5+ic^ + ^^. (5.13) 



Here we have used Lemma 14.11 Corollary 15.11 the first condition in ( 15 .6 aft and Cauchy esti- 
mates. We then estimate p + . From (I5.12p and ( 15.5b|) we obtain 

\\p+\\v-£ < 4 d z H \\v-Z/2,zM\ d w(\\v-tM\Po\\v-z/2 (l + e||^C||^/2||/3i|U-?/2) 
+ e\\d w Ho( 1 \l_s\\p 1 \\ v _s /2 

_ 2 2K5 AC D K5 ( _ 1 2K5 AC D K5 \ _ x . AC D K5 
<e2k c h—z — (l + e« — — J+€K \\d w H o Ci\\v-t — 7 — 

< AeC D K (2k~ 2 C h (l + 2eC D K~ 1 ) + k-^H o Ci||^)r^ 

Now, introduce 

c§ = AC D K (^2k- 2 C h (1 + 2eC D K^) + kT 1 (c' h + 6+ ~C A K6 + ^^H ) 
so that (I5.6ep holds true. 
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For A + we have using (15.6b|) : 

A + (w + ) = d 2 z+ H + (w + , 0) = d 2 z H{d z+ zf + d 2 w H(d z+ w) 2 + 2d z d w H(d z+ z)(d z+ w) 
+ d z H8 2 z+ z + d w Hd 2 z+ w 

= d 2 z H(I + edvtxhf + £ 2 dlHf3 2 + 2ed z d w H(I + ed^ih)^ 
+ ed z H{d 2 w Uh) 2 + d w (A(w + , 0) + ed w H/3 2 (w + , 0). 

Here we used (I5.3al) . (15.3b[) and (I5.3cp . and all derivatives of H on the right hand sides are 
evaluated at (w,z)(w+,0) = (w+,Ci), cf fOal) . flO) . ( Ojl . Using d 2 H(w + Xi) = A (™+) + 
d 2 z r(w + X\) gives, with (15. 7p . 

— A(w + ) = e / d w A(w + + esf3 )(3 ds + d 2 r 
Jo 

+ 2ed 2 z H(d w (A) + e 2 d 2 H(d w (A) 2 

+ e 2 ^^ 2 + 2ed z d w H(I + ed w CMPi 

+ e«9 z #(<9 2 G(A) 2 + d w (ik) + e0 w iffr. (5.14) 

Since r is cubic in z by (H3) we have d 2 r(w, z) = \ (1 — s) 2 d^r(w, sz)zds. Therefore, apply- 
ing Cauchy-estimates on d w A and d z r and using Lemma T5.14 (I5.5a[) . (I5.5bj) . and Lemma I4TT1 
we obtain 



2C A AC D K5 6C r 
\\A+ - A\\ u -s < e— + — Kb 



2C H ( 2K5AC D K8 , f2K5AC D K5 x2 



k 2 \ £ «£ 



( 2K8AC D K8 \ 

\1 tfT) 



2 2||^FoCi||^ f4C D K8\ 2 2 .. .. / 2AT54C D A:5\ 4^X5 



+ 4 (^) 2 + ^^) +£ ||^o Cl ||^, (5,5) 

so that, by (15.131) and conditions (15.6a[) . 

P + -A||^<c Ail 5 + ^^ (5.16) 



for e sufficiently small. Here 

ca,i 

and 



5 : 



32C // C , z)K 2 <5 32eC fc Ci,AT J j. 32C H C 2 D K 16C H C D K 2 5 



CA,2 = SC A C D K + 3 + I * + 3 + 

flQeClK 8C D Ku ( 2C D \ 16C D Ku \ 
+ — -t— + 5 — 1 + + 3 \\d w H o Cil^-e, 
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and we use that £ < v and £ > 2Kb. The rows in H5 . 1 5 j) correspond to the rows of A + — A 
as they appear in (I5.14p . As a consequence we obtain 



^\\y\\<\\A + \\\\y\\ + <^ + c A , 2 )5\\ y \\, 



and so, if we choose e > small enough such that 
bearing in mind that 5 = 0(e) by (H2), then 



(5.17) 



k\ 2 v e 



^--(^f + c A , 2 )5)\\y\\<\\A + \\\\yl 



so that (15.6dj) holds true. Now finally for r + , using (15.6bj) . (15 .ip and (I5.3aj) we obtain 

r + (w + , z+) = H(w, z) - (h + (w+) + (p+(w+), z+) + ^(A + (w + )z + , z+) 
= h(w) — h + (w + ) 

+ (pH - P+(w+), z+) + (p(w), C(w)) 
+ ^((A(w)-A + (w + ))z + ,z + ) 

+ (A(w)C(w), z + + -C(w)) + r(w, z), (5.18) 

Here we used that 

~(Az + ,z+)-~(Az,z) = (A(z+-z),-(z + + z)) 
We now write h(w) — h + (w+) as (h(w) — h(w + )) + (h(w + ) — h + (w + )), and, using ( I5.3aj) . write 

h(w) — h(w + ) = h(w + + e/3) — h(w + ) = / d w h(w+ + se(5)ef3ds. 

Jo 

Similarly, we split A(w) — A + (w + ) and r(w, z) — r + (w + , z + ), writing 

A(w) - A(w + ) = A(w + + e/3) - A(w + ) = e f d w A(w + + se(3)f3ds, 

Jo 

and 

r(u>, z) — r(w + , z + ) = r(w, z) — r(w, z + ) + r(w, z + ) — r(w + , z + ). 
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We then use (15.6cj) to estimate A{w + )— A + {w + ). Using (15. lip . (15.5a|) . Lemma ISTT] Lemma |47T| 
and (I5.6ep we obtain 

e5 „,AC D K5 
\\r+ - r|| v _ €i<T _£ < c h>1 6 + c hy2 -^ + eC h 

+ C D (5 + 5 + ) + K5 2 

1 2C A AC D KS 2 1 fec A>1 6 \ 2 



l r/n C r AT<5 8C r C D KS 
2 K5) + — + e — 



+ C A K8{C D + -K8) + + e^T^ (5.19) 



< c ri i<5 



where 



Cr,i = c M + C D + #5 + -c^Cf, + CUC D A + \c A K 2 5 + 

c r , 2 = c M + 4C D KC' h u + C D c^ + AC 3 D C A K + ^c A)1 C£ + 80-^. 

As above, the rows in the estimate of r + — r in (15. 19[) correspond to the rows of r + as they 
appear in (I5.18p . 

Finally, we note that ( 15.1 ip gives 

\\h+\\„-s < C h+ = C h + c h>l 8 + °h,2-^ 

and 

||9«,/i + || v -3e/2 < \\d w {h+ - /i)|U-3f/2 + \\d w h\\ v ^/ 2 

using a Cauchy estimate. ■ 

Next, we assume that (H4) holds true. The following lemma shows that the equilibrium 
is contained within the improved slow manifold {z + = 0}. 

Lemma 5.4 Assume (H4)- In the (w + , z + )- coordinates the equilibrium is given as (u>+, 2+) = 
(w e ,0). 

Proof First note that z e = implies that C{w e ) = 0. By the definition of the symplectic 
transformation from the Iterative Lemma 15.31 and the local uniqueness of the equilib- 
rium, assumed in (H4), it will suffice to show that (w+, z+) = (w e , 0) satisfies z\) = 
(w+,z+), i.e. solves the following system: 

4 = 2 e -C(« e ,0 = -CK,<), 

u\ = u e + e{d v+ Cy, x e ) - e(d v+ ( x , y%) = u e - e(d v+ ( x , y e + ), 
v e + = v e - e{d u ( y , x e ) + e(d u C x , y%) = v e + e(d u ( x , y e + ). 

Insertion shows this directly. ■ 
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To prove Theorem 13.21 we proceed as in the proof of Theorem l3.lt We will first apply the 
transformation: 

: (V + i(v - Co)) x (S + i(a - &)) (V + i(u - £ /2)) x (S + z(a - £o/2)), 

(u>l,Zl) (lUo^o), 

once with £ = £o = 0(1) < max (j/o — u, e>"o — <I_) large. By the Iterative Lemma [5.3[ using 
that 1 1 Jo 1 1 = 0(e) by (H2), we can therefore take 8% = C$e 2 , as in (I4.16p . and start from 
(V + zVi) x (S + iai) with v x = ^o-£o an d <J\ = fo-£o- Also, ||/i x - h \\ ui , \\A 1 — AqW^, ||r x - 
r o||yi,tri = ^(e)- We then continue to apply the transformations 

^+ :(V + iv n+ i) x (5 + za n+1 ) -)• (V + z(v n - £ n /2)) x (5 + z(a n - f n /2)), 

(u>„ + l, Z n+ l) H> (!D W Z,,), 

successively with £ n = 2emax{c,5 n , i^ n } for n > 2 and £ x = 2emax{cj 1 , iTi, y/2Ki CdC$}, so 
that Lemma [5.21 applies and 

<Wl < ^n, (5.20) 

for 1 < n < N. Here 5 n = \\p n \\u n with v n = v x - Y^lZ\ih o n = o x - YTiZ\ii- Then 
<Wi < 2 _n <5 1 . To control the possible growth of the constants Ch n ,C' hn ,CA„,C rn , K n with 
n we may first restrict N not only by the conditions that the domains remain non-empty: 
v n — 6v > E_ and <7jv — £n > <L but also by 

C^ n , C An , C r „,^ < 2max{C f fcl , C^, C Al , C ri , K,} for all n < N. (5.21) 

Then we can estimate the constants Ch,i, Ch t 2, c$, ca,i, c a ,2 and c rj i, c rj 2 in the Iterative 
Lemma 15.31 uniformly with respect to n < N . This allow us to estimate the growth of 
the constants Cfr n ,C' h ,C An ,C Tn , (summing over n, using 5 n = 2~ n+1 6\, 5i = 0(e 2 ) and the 
geometric series formula, similarly as for the estimation of f3 n in the proof of Theorem 13. 1[ 
see (I4.26p . (I4.27P ). Moreover we estimate K n using ( 15 .6 dp . similarly as we did for a n in 
( I4.22p . (I4.23p . From this we are directly led to the conclusion that the conditions ( I5.2ip pose 
no restrictions on N, as they can be satisfied uniformly in n by choosing e sufficiently small. 
Furthermore, £ n < £ for some £ = (9(e) giving N = 0(e~ 1 ) possible step. Since 5\ = 0(e 2 ) 
we finally conclude that there exist new constants so that 

\\h N+1 - hi\\„ < c h e 2 , \\p N+1 \\u < c 8 e 2 e~ c/e , 
\\A N+ i - < c A e, \\r N +i - nW^a < c r e, 

for e < e. This completes the proof of Theorem 13.21 

6 An invariant two-dimensional slow manifold 

A consequence of Theorem 13.21 is the persistence of a two dimensional normally elliptic slow 
manifold with exponentially small gaps. We will sketch this result here which extends a 
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result in [9] to several fast variables. Consider a real analytic slow-fast Hamiltonian system 
H = H(wq, Zo) with 1 slow degree of freedom, dz fast degrees of freedom and assume that 
zq = is invariant for e = for w G V C M 2 . We assume that {zq = 0, e = 0} is filled with a 
non-degenerate family of periodic orbits. Then by Theorem 13. 21 there exists a transformation 
(w, z) H- (u> , Zq) that transforms H system into 

H = h{w) + (p(w), z) + -{A(w)z, z) + r(w, z), 

with p = 0(e~ c / £ ) provided e is sufficiently small. The Hamiltonian system h = h(w) is inte- 
grable since it is a one-degree of freedom system. By assumption we can therefore introduce 
action angle variables (<f), I) w so that h = h(I) with djh ^ 0. This transformation does 
not depend upon the fast variables and can therefore directly be lifted to the full space to 
give: 

H = h(I) + (p(0, /),*) + ^(A(cf>, I)z, z) + f(cf>, I, z), 

where p(0, 1) = p(w), A((f>, I) = A(w) and f(0, /, z) = r(w, z). We will henceforth suppress 
the dependency on I. Since dj h ^ we can solve the equation H = E of energy conservation 
for / = Ke((/),z) when z is small. Differentiating H(<f>, Ke, z) = E with respect to z gives 
d z K E = —(diH)~ l d z H and therefore: 

ed^z = -J z d z K E , 

where J z is the matrix from ( 13. 51) . Here 

K E = h~ l (E) - (p E , z) - ^(A E z, z) - r E , 

with p E = {d I h)- 1 p = £>(e- c / £ ), A E = {d I h)~ 1 {A + 0{e~ c ' e )) and r E = 0(z 3 ). Fix e 
small and introduce p? = supp^ = 0(e~ c l e ) so that pe = p 2 p with supp = 1. Also let 
p 2 B = A E — A, i/j = (<9//i) -1 and z = pz. We then introduce 

Q = Q(if>, z) = /i(p, z) + ((A + p 2 B)z, z) + pf, 

f((f),pz) = p 3 r E (ipyZ) so that 

ed^z = J z d z Q = pJ z (p + pBz + d z f) + J z Az. 

Let 

P, : = 0} -> = T}, T = 2 7 r(9 / / i )- 1 , 

be the corresponding stroboscopic mapping. It is symplectic since the system is Hamiltonian. 
Also .Po(O) = 0. The persistence of this fixed point for p ^ provides the persistence of 
the periodic orbits, which we have parametrized by E. For this we consider the linear 
map 8 z Pq = ^(T), the monodromy matrix, where ^ = ^(ip) satisfying ^(0) = /, is the 
fundamental matrix of the system 

ed^z = J z Az. 
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The eigenvalues of \I/(T), Xi(E), . . . , \2d z (E), are the characteristic multipliers and they 
depend upon E. We can apply the implicit function theorem to the set 

£ = {E\\\{E) - 1| > y/fi,i = 1, . . . , 2d z }, (6.1) 

since the perturbation is O(fj) = 0(er c ^ 2e '). Generically one would expect the singular 
set {E\\\i(E) = 1 for some?} to be discrete containing simple solutions and without accu- 
mulation points. In the affirmative case the complement of the set S from ( 16. II) would be 
exponentially small, and the slow manifold {^o = 0, e = 0} would therefore perturb to an 
invariant slow manifold {zq — e < eo with exponentially small gaps. Note that this 

result does not depend upon resonances between the fast variables. 

7 Conclusion 

We proved the existence of slow manifolds which are invariant up to an exponentially small 
error in general and Hamiltonian analytic slow-fast systems. The slow manifolds were con- 
structed so that they included nearby equilibria. The approach we used for general systems 
is due to R. S. MacKay and does not require the transformation of the slow variables, and as 
a consequence, our results in this case hold true for unbounded fast vector-fields. In future 
work we aim to extend our results to unbounded slow vector-fields and infinite dimensional 
slow-fast Hamiltonian systems. Here we believe that there is great potential in combining 
the methods used here with the methods due to K. Matthies and A. Scheel [T9] . 
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